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ABSTRACT 


In situ measurements of salinity and temperature fluctu- 
ations at depths to 14 meters indicate distinct dependences 
at different times of the day. The variance of the salinity 
fluctuations decreased with increasing depth, but was greater 
just after sunrise than just prior to sunset. The variance 
of the temperature. fluctuations decreased with increasing 
depth just prior to sunset, but increased with depth immedi- 
ately after sunrise. The correlation length of the sound 
index of refraction was calculated by using the variance of 
the sound velocity fluctuations, and the variance of sound 
amplitude modulation in the theory of Mintzer. This analysis 
shows that microstructure patch size increases approximately 
linearly with depth. The power spectral densities of the 
salinity, temperature and sound velocity fluctuations show 


peaks of energy corresponding to dominant ocean wave 


frequencies. 
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I. INTRODUCTION 


A. HISTORY 

The study of sound propagation in the ocean and the ef- 
fects of oceanographic parameters (particularly temperature) 
on it has been the subject of considerable research since 
early World War II. Wilson [Ref. 1] formulated an equation 
describing the relationships between sound velocity, tem- 
perature and pressure in distilled water, and then extended 
his work to an oceanic regime by utilizing filtered seawater 
diluted with distilled water to obtain desired salinity 
dependence [Ref. 2]. Since Wilson's efforts in this field, 
many investigators have attempted to revise the sound speed 
equation with newer techniques and theory. The reasons for 
revised formulations included, among others, the use of new 
type velocimeters [Ref. 3] the suspicion that standard sea 
water did not maintain a standard salinity to sound speed 
relationship [Ref. 4], the desire to obtain simpler expres- 
sions, and the use of a more realistic range of values for 
the oceanographic parameters [Ref. 5]. 

While these investigations into the parametrical rela- 
tionships to sound velocity were proceeding, other investi- 
gators attempted to measure the magnitude of the fluctuations 
of the important parameters affecting sound propagation in 
seawater.  Liebermann [Ref. 6] studied the effects of tem- 
perature inhomogeneities on sound propagation. He found that 


thermal inhomogeneities do have an effect on the scattering, 


Nen 





intensity, and refraction of the sound transmission. Shon- 
ting [Ref. 7] also measured the thermal microstructure, his 
results indicating the importance and magnitudes of tempera- 


ture fluctuations in the ocean environment. 


B. SIGNIFICANCE OF PROBLEM 

The characteristics of sound propagation in a deep ocean 
environment (defined in this paper as the depth below the 
thermocline) are well known and predictable. This is due 
to the fact that in the deep water regime the fluctuations 
of temperature and salinity are minor. As a result of this 
predictability, sound is used properly and constructively 
in deep waters. For the United States Navy this specifically 
relates to optimum usage of sound in the field of anti- 
submarine warfare. 

However the propagation of sound in the shallow ocean 
environment and in the near surface region is not as pre- 
dictable, and as a result seriously hinders the user and 
the planner in their efforts to utilize sound to its best 
advantage. The near surface region usually plays a part in 
long range propagation of sound. 

Sound velocity amplitude and phase in the region between 
sound source and sound receiver change in the near surface 
regions. They change due to small scale variations in oceano- 
graphic parameters including possibly the presence of micro- 


bubbles and organic content in the regime. 
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e. OBJECTIVES OF THIS RESEARCH 

The main motivation in this research was to statistically 
describe the oceanic parameters of temperature and salinity, 
and to relate their statistics to the statistics of sound 
velocity in a near surface environment. It was hoped to ob- 
tain enough data to enable one to draw conclusions concern- 
ing the variation of the fluctuations with respect to depth 
and time of day. If such concrete conclusions could be 
reached in the near surface region, then the users of sound 
either in oceanographic research or in present anti-submarine 
warfare would be afforded new tools with which to update the 


state-of-the-art in sound use and understanding. 
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II. EXPERIMENTAL PROCEDURE 


A. OCEANOGRAPHIC ENVIRONMENT 
I Location 

Due to its E pi location of being rigidly 
fixed in a shallow water environment the N.U.C. Oceano- 
graphic Research Tower, approximately one mile off Mission 
Beach, California, was chosen as the location for the pro- 
ject. It allows for continuous data accumulation and 
monitoring as well as some analysis in a near laboratory 
controlled environment. 

Figure 1 is an illustration of the tower. It is 
fixed by supporting pins driven 63 feet into the ocean 
floor. Electrical power is supplied from shore thereby 
guaranteeing stable voltage and frequency. The tower is 
located in 60 feet (18 m) of water with a sandy bottom 
and is free of water traffic most of the time. The upper 
level of the two level tower was the location for elec- 
tronic support equipment. The lower was for handling 
equipment. 

2. Time of Year 

The experiment was conducted from the afternoon of 

me October 1971 through the morning of 22 October 1971. 
3. Meteorological, Oceanographic, and Biological 
"Conditions at Time of Project 
a. Weather State: clear to light haze 
b. Wind: 21 October - 300-320 "From NEO T R ES 


22 October - 065-110 from 2 to 8 kts 
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NUC Oceanographic Research Tower. 





c. Wind Waves: 0 to % ft from crest to trough 
d. Swell: 2 ft from Northwest 
e. Bubbles: 1 square foot in vicinity of tower 
stanchions during all runs. During run #4 bubbles were 
present on the surface in patches of size 20 ft by 2 ft 
Hear vicinity of tower. 
f. Biologics 
(1) Seaweed. Small amount of sea grass on 
surface during run #4. 
(2) Animal Life. Seal in vicinity during run 
#6, porpoise and seals in vicinity during runs 7 and 8. 
g. Temperature Structure: thermocline varied from 
depths between 40 ft to 60 ft below the surface (see 
Big. 2) 
h. Salinity: Nansen Cast results, as determined 


from Hytech model 621 salinometer, found in Table III, 


B. EQUIPMENT USED DURING EXPERIMENT 
l. Bissett-Berman STD Model 9006 
a. Temperature Sensor 
(1) Time Constant: 0.35 see 
(2) Output: 0 =- 10 mv DC 
(3) Temperature Range Used: 14 - 19° C 
(4) Accuracy: :0.02° C 
b. Salinity Sensor 
(1) Time Constant: 0.35 sec 


(2) Output: TOFS OE DC 
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Bathythermograph Traces. 
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(3) Salinity Range Used: 37.5 — 39.5 ppt 
(before correction) 
(4) Accuracy: +£0.03 ppt 
c. Detailed Description Found in Appendix A 
2. Thermistors (3) 
a. Time Constant: 150 milliseconds 
LD; Output: 060553 Les perc 
c. Temperature Range: Thermistors nulled using 
Wheatstone Bridge Circuit prior to start of each run. 
d. For further description see NPS thesis of LCDR 
Duchock, March 1972. 
3, Ramsay Corporation Mark I SVID 
a. Time Constant: 160 microseconds 
bD. Output: ONO y DE 
c. Sound Velocity Range: 1400 - 1600 m/sec 
d. Accuracy: +.01 m/sec 
e. For further description see NPS thesis of LCDR 
Duchock, March 1972. 
4. Wave Height Sensors 
a. Baylor Wave Profile Recorder System 
(1) Output: 50 millivolts per foot of height 
(2) Accuracy: 1% or reading 
(3) For further description see NPS thesis of 
mee Bordy, March 1972. 
b. Interstate Electronics Corporation Pressure Wave 
Gauge 


(1) Output: .25 v/psi 
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(2) Accuracy: 2% of reading 
3.  Parerneke Velocity Sensor 
a. Engineering Physics Company Electromagnetic 
Current Meter 
(1) Output: variable voltage per meter/second 
(2) Accuracy: % of full scale reading 
b. For further description see NPS thesis of LT 
Bordy, March 1972. 
6. High Pass Filter #1 - Krohn-Hite Model 3340 
a. Gain: OdB 
b. High Pass Filter Cutoff: 0.01 HZ 
7. Band Pass Filter #2 - Krohn-Hite Model 330-A 
a. Gain:  OdB 
b. low Pass Cutoff: 0.02 HZ 
High Pass Cutoff: 2000 HZ 
8. Band Pass Filter £3 - Velocimeter 
a. Gain: OdB 
b. Low Pass Cutoff: 0.02 HZ 
High Pass Cutoff: 2000 HZ 
9. Amplifier #1: Preston Scientific Model 8300 xwb 
a. Gain: 500 
10. Am 


lifier #2: Preston Scientific Model 8300 xwb 











a. Gain: 1 

ll. Amplifier #3: ' Hewlett-Packard Model 2470 
a. Gain: 30 

12. Amplifier #4: Hewlett-Packard Model 2470 


a. Gain: 100 


1.9 





153; Amplifier #5: Hewlett-Packard Model 2470 
a. Gain: 100 


14. Amplifier #6: ' Hewlett-Packard Model 2470 
a. Gain: varied 100-300 
15. Sangamo Magnetic Tape Recorder Model 3500 
a. :14 channel with 11 channels allocated for 
oceanographic parameter data were utilized. 
b. Channel Allocation: 
(1) Ramsay Temperature Sensor (not used) 
(2) Ramsay Velocimeter 
(3  Bissett-Berman Salinometer 
(4) Turbulent Velocity (horizontal component) 
(5) Turbulent Velocity (vertical component) 
(6) Baylor Wave Height 
(7) Sound Amplitude Modulation 
(8) Pressure Wave Height 
(9) Thermistor #1 
(10) Bissett-Berman Temperature (not used) 
(11) Thermistor #2 
(12) Thermistor #3 
(13) Sound Phase Modulation 


i (14) Voice (not used) 


C. EXPERIMENTAL CONFIGURATION 
l. Schematics 
a. Bissett-Berman STD - See Fig. 3 
b. Ramsay Velocimeter - See Fig. 4 


C. ‘Temperature Thermistors - See Fig. 5 
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FIGURE 3. Bissett-Berman STD Flow Diagram. 
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FIGURE 4. Ramsay Sound Velocimeter Flow Diagram. 
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FIGURE 5. Thermistor Flow Diagram. 
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2, Geometry of Sensors 


The sensors' geometric relationships were deemed to 
be of major importance in view of potential acoustical and 
physical interferences which could have had an adverse ef- 
fect on the experimental results. In addition to the attempt 
to describe the oceanic regime in terms of the oceanic 
parameters two acousticians simultaneously conducted ex- 
periments related to sound amplitude modulation, sound 
dispersion, and phase fluctuations. (See NPS thesis of 
LCDR W. J. Smith, Jr., December 1971 and thesis of LCDR 
Juergen Rautmann, December 1971.) 

In order to reduce the effects of mutual interfer- 
ence between the sensors especially in the path of sound 
propagation between hydrophone and receiver, a physical 
separation of the sensors was mandatory. The geometrical 
separations between the sensors as shown in Figs. 6, 7 were 
used. Photographs of the Bissett-Berman STD and the geo- 


metrical relationships of the other sensors to it appear 
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III. DATA COLLECTION AND RUN DESCRIPTION 


Data collection and run description will be found on page 56 under 


Appendix B. 


29 





IV. DATA REDUCTION 


A. INTRODUCTION 

The analog data were recorded for each run on the 
Sangamo Model 3500 magnetic tape recorder operated at a 
speed of 1 7/8 ips using FM electronics. The recorded 
data were then transferred to a Brush Mark 200 Strip-Chart 
recorder. In some instances, in order to produce better 
resolution, the recorded data were transferred to a Brush 
Mark 2 Strip-Chart recorder through a Krohn-Hite Model 
3340 filter set on low pass-max flat (see Fig. 12). The 
printed analog data were converted to digital data on 
magnetic tape using the Fleet Numerical Weather Central, 
Point Pinos, California tracing digitizer. The taped data 
were transferred to punched IBM data cards compatible for 


statistical analysis on an IBM 360-67 digital computer. 


B. MAGNETIC TAPE TO STRIP-CHART RECORDING 

The recorder playback speed was seven and one half 
inches per second, which is an increase of four times over 
the data collection speed. The strip-chart recorder was 
run at a speed of five millimeters per second. The speed 
differentials between the record, playback, and strip-chart 
phases were programmed in the IBM program used for data 
analysis. The sensitivity setting on the strip-chart re- 
corder was important for data presentation and analysis. 


The sensitivity settings varied from run to run and with 
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FIGURE 12. Representative Oceanographic Fluctuations. 
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each parameter. The settings used can be found in 


Appendix B. 


C. ANALOG TO DIGITAL TAPE 

In order to utilize the IBM 360 computer to determine 
the statistics of the fluctuations, the analog data on et 
strip-chart record had to be converted to digital form. In 
order to effect the conversion use was made of the Calma 
tracing digitizer at Fleet Numerical Weather Central's 
Bunt Pinos, California facility. The tracing digitizer 
transfers analog data to digital data on magnetic tape at 
a sampling rate of 100 samples per inch in both the x and 
y directions. The length of each run analyzed amounted to 
approximately the first nineteen minutes of the run while 
the number of data points varied depending upon the posi- 
tioning of the strip chart on the digitizer. It was im- 
posant for future correlations that the tracing of all 
parameters began at the same starting point. This was 
ensured by starting at the sharp line transient signaling 


G Start of a particular run. 


D. CONVERSION OF TAPED DIGITAL DATA TO DATA CARDS 
Conversion of taped digital data to a form compatible 
for use on the IBM 360 Digital Computer, i.e. data card 
form, was accomplished using Fleet Numerical Weather 
Central's Computer Facility. Fourteen data points were 


recorded on each IBM card. The number of data points per 
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run varied depending on the positioning of the strip chart 
on the Calma Plotter at Point Pinos, but the average run 


produced approximately 5400 data points. 
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V. DATA ANALYSIS 


PON INTRODUCTION 

Attention was focused entirely on analysis of the 
statistics of the fluctuating parameters. Two arguments 
led to this method of analysis: first of all, the statis- 
tics of the oceanographic variables are important in deter- 
mining the statistics of sound propagation. Since the 
major objective of the experiment was to define the medium 
in order to determine its effect upon the propagation of 
sound the statistical method of analysis was the natural 
approach. The second reason for using statistics as the 
major tool for analysis was due to the relatively wide 
separation between sensors which was required to prevent 
instrumental interactions. In order to enable one to dis- 
cuss the interrelationships between spaced sensors it was 
necessary either to correct for the lack of idential point 
location or to consider the statistical relations, assum- 
ing that the medium was statistically homogeneous in the 


region of the sensors. The latter is the easier choice. 


B. STATISTICAL PROGRAMMING 


The programs used in the data analysis appear in 


Appendix C. 

Program A analyzed the data for mean, variance, normal- 
ized autocorrelation function, and power spectral density. 
The standard Blackman-Tukey method was used to compute 


power spectral density with a Parzen window as the spectral 
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window. Program A was modified to compute variance by two 
methods. The first method calculated variance by integra- 
ing the power spectral density. A subroutine (aver) was 
then added to calculate variance directly from the time 
series data. 
Program B was used to output data from which a histogram 

could be plotted of the number of times a particular value 
of a parameter occurred versus the value of that parameter. 


In utilizing program A several important inputs were 


necessary: 

NTS: number of data points 

MLAG: number of time lags = 10% (NTS) 

DT: time increment in seconds/sample. This parameter 
took into account the different driving speeds of the strip- 
chart recorder, playback and record modes of the tape re- 


corder, and the analog-to-digital digitizer. 


1 
DT = 
inch 100 samples 
(5/4, SSA A a 
25.4mm neh 


0.2032 sec/sample 


T 
Af: frequency increment 


2 MLAG (DT) 
FBHZ: lowest frequency 


FEHZ: highest frequency = Nyquist Frequency 


1 
EE A TZ 
2(.2032) 


Calibration factor: variable on each run 
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The raw data on the IBM cards had the units o£ inches. 
The calibration factor converted inches to volts per milli- 
meter by incorporating the Brush recorder sensitivity set- 
ting into the program input thereby giving the data units 
of volts per inch. When each data point was multiplled by 
the calibration factor the results produced were in voltage 
units. 

| $ lee 25.4 millimeters volts 

SEE a EN 

millimeter inch inch 

C. STATISTICAL ANALYSIS 

m” Introduction 

In order to convert the program output in volts to 

the correct oceanographic units and values, conversion 
factors for sensor voltage sensitivities, amplification 
factors and filter gains were used. On location monitoring 
of the data combined with post-experimental review in- 
dicated some runs higher in quality and accuracy than 
others. For such reasons runs 3,4,5,6,7,9,10 were analyzed 
for salinity and temperature, and runs 4,5,6,8 for sound 
velocity. The conversion factors used were: 

Lthermistors: 5.53 v/*C 

Salinity: 2.5v/% 


m 
sound velocity: 1.5 v/— 


sec 
2. Variances of Oceanographic Fluctuations 
The values of the variances of the oceanographic 


parameters as computed by integrating the power spectral 


densities are indicated in Table I. 
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The values of the variances of the oceanographic 
parameters as computed directly from the time series data 
are listed in Table II. 

It was noted that there was only a very slight dif- 
ference between the values of variance computed by both 
methods. However the variances used in subsequent analysis 
were those computed directly from the data thereby ensuring 
a higher degree of accuracy. Additionally ‘all subsequent 
analysis’was performed using temperature data measured by 
thermistor #1 since the quality of the fluctuating signal 
from that sensor was the best and most consistent throughout 
all runs. (Thermistor #2 data used during run 4 since 
thermistor 2 was operating reliably at that time). Figures 
13-15 are graphical displays of variances versus time of day 
with depths at which runs were made noted in parentheses. 

A survey of the literature indicates that this ex- 
periment was unique in its attempt to establish statistical 
relationships between the described oceanographic parameters. 
As a result little direct comparison with previous results 
is possible. 

3. Variance of Temperature Fluctuations 

Two past investigators have studied temperature fluc- 
tuations in the near surface region on approximately the same 
time scale as was used in this experiment.  Shonting used a 
double-beaded thermistor manufactured by Fenwall Electronics, 
Inc., Framingham, Massachusetts, to observe thermal micro- 


structure. He conducted his measurements in the Bahamas in 
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Sound Vel. | Temp #1 Temp #2 Salinity 
Run |Time |Depth| Variance Variance | Variance Variance 
21 OCT 
3 J1616 | 9.0m SEET 
PODES 


1648 l13.9m| .0294 (M/S) ?* Í 
.0003 .0230 . 0002 .0023 | 


m „0072 “0052 
om -0023 „0013 


10__|0802 m 


TABLE I 
z 00 ji 6 (%o 


E LL. Sound Vel. | Temp #1 Temp #2 
Run |Time |Depth| Variance Variance |Variance 
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3. |1616 | 9.0mj - 0143 (c0 eelere eC) * 
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4 NT 13.9m| .0293(M/S)*| .0129 = .0012 
5 n 6.9m| .0003 |  .0184 . 0002 .0023 
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August 1967. Analysis of his measurements taken between 
1323-1400 indicated representative values of variance of 
approximately 25 x 107°(°C)? at two meters depth and ap- 
proximately 72 x 107°(°C)* at 20 meters, a decrease of two 
orders of magnitude with depth down to 20 meters. Although 
the variance was much greater in the present experiment 
than in Shonting's, there was in the late afternoon period 
a one order of magnitude decrease in temperature variance 
with depth down to approximately 15 meters. 

Liebermann also recorded temperature fluctuations in 
the near surface region at a depth of 50 meters. A rough 
calculation of variance from Liebermann's Fig. 1, which was 
given as a representative temperature record for July and 
August (time of day unreported) for the continental waters 
from Southern California to Alaska, resulted in a variance 
of approximately 17 x ONE This figure is somewhat 
closer to data of Fig. 13 but the lack of time of day in- 
formation and depth dependence prevents conclusions from 
being drawn. 

Due to the maximum temperature gradient being effected 
just prior to sunset the strongest fluctuations of tempera- 
ture were noticed at all depths at that time. The variance 
decreased with increasing depth. However in the short 
period just after sunrise the temperature fluctuations were 
noticed to increase with depth. Sagar [Ref. 8] quotes, from 


a personal communication, Deacon, of the National Institute 


a \ = un S - E. d e 
È u u A B Š £ q M xu cn = TW LA mm 
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FIGURE 14. 


Plot of Salinity 
Variance Versus Time 
of Day (Depth in 
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FIGURE 15. 


Plot of Sound Velocity 
Variance Versus Time 
of Day (Depth in 
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This behavior could possibly have been caused by the "con- 
vective lag" phenomenon. James [Ref. 9] describes this 
behavior as occurring in the early morning just after sun- 
rise. Insolation increases to the point where it is equal 
to the outgoing heat loss from the water. However heat 
loss occurs exclusively at the surface, with consequential 
continued cooling of the surface. As insolation increases 
in the hours just after sunrise, surface cooling continues 
while the lower levels are slowly heated. Eventually, of 
course, the lag is overcome and surface heating begins 
causing growth of small cells with large fluctuations near 
the surface and large cells of small fluctuations at greater 
depths, as described by Deacon. 

4. Variance of Salinity Fluctuations 

The analysis of the plots of salinity variance versus 

time showed a decrease in salinity fluctuations with depth 
during the late afternoon period. The fluctuations were 
smaller during the period just prior to sunset than just 
after sunrise (assuming that the 7.8m reading is in error). 
Wi-Wwdccrease in variance of salinity with increasing depth 
in the late afternoon was expected ee a stronger salinity 
gradient would exist in the surface layer where evaporation 
was initiated. However, it is speculated that convective 
mixing caused a more nearly isohaline condition to exist in 
the water column by sunset than existed just after sunrise 
thereby accounting for the difference in magnitudes of 


variances at those two times of the day. 
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In the section "Parametric Interrelations Between 
Salinity, Temperature, and Sound Velocity" the theoretical 
calculation of sound velocity variance is discussed, and 
further analysis of the difference between theoretical and 
experimental results is made. 

5. Probability Density Function Results 

Salinity and sound velocity information were put into 
Program B in order to establish the type of distribution of 
the populations sampled. In all cases the distributions were 
Gaussian in nature. Figures 16 and 17 are examples of the 
distributions. 

6. Normalized Autocorrelation Function 

The normalized autocorrelation functions of the oceano- 
graphic parameters (see Appendix D) indicated a long correla- 
En time for temperature data thereby indicating relatively 
well mixed thermal structure. Conversely the relatively low 
degree of correlation for both salinity and sound velocity 
indicated an inhomogeneous medium with respect to those two 
parameters. The peaks and troughs found on the salinity and 
sound velocity autocorrelation plots indicated the presence 
of patches of higher and lower salinity water as well as 
patchy areas of water yielding higher and lower values of 
sound velocity. 

ie Power EE EE 
Plots were made of power spectral densities of salin- 


ity, temperature, and sound velocity fluctuations (see 


Appendix E). 
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Salinity Histogram. 


o 
r 
m 
O 
H 
E 
A _ nn 
O O © 
LO) O LO) 
N N ri 


Ue G$ (¿ATO?) UT SATRA SIUETIEA EE EES EK 


45 


O 
O 
= 
ed 


O 
LO 


JO S8DUEAANDDO JO £d 


2509 Feoi or perce ma 


"BU 


Salinity Variance in (volts)? 





Sound Velocity Histogram. 
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The nature of the plots indicated peaks of "energy" 
centered about frequencies of 0.17 and 0.08 Hertz. This 
coincided with the dominant wave frequencies existing at 


the time of the experiment. 
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VI.  PARAMETRIC INTERRELATIONS BETWEEN SALINITY, TEMPERATURE, 


AND SOUND VELOCITY 


A. INTRODUCTION 

With hopes of shedding light on possible heretofore un- 
known variables capable of affecting sound velocity a com- 
parison was made between sound velocity variance derived 
from Wilson's sound speed equation and that measured experi- 


mentally with the Ramsay Velocimeter. 


B. DERIVATION OF THEORETICAL SOUND VELOCITY VARIANCE FROM 
WILSON'S EQUATION 


1. Wilson's Equation 


C (t) = 1449.2 + 4.623 T (t) - .0546 T? (t) 
+ 1.391 (S (t) - 35) 


2. Computation of Variance 


Nat C = 0 + (4.623) var £ =-(.0546) kar T) 
+ (1 391), 2 len 
3. Dropping the Term, var (T*), as Negligible Yields 
the Final Result 


var C = 21.372 var T + 1.935 var S 


The utilization of the above formula yielded the 
following results, with the velocimeter results shown for 


comparison: 
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RUN VARIANCE Ctheory VARIANCE o o 
3 .3095 (m/sec)? - 

4 .2770 (m/sec)? .0293 (m/sec) ” 

5 .3979 (m/sec)? .0003 (m/sec) ” 

6 .0968 (m/sec) * .0072 (m/sec)? 

g 0437 (m/sec)? .0021 (m/sec)? 

9 NS E a E 
10 0273 (m/sec). S 


It was hoped that the results would indicate a theo- 
retical variance slightly smaller than experimental variance 
thereby implying the presence of other parameters such as 
bubble presence to which the 3 MHz acoustic velocimeter 
would have been insensitive. However, as seen above, the 
experimental variance was much smaller than that calculated 
from the variances of temperature and salinity. Analysis 
of the wide discrepency led to the conclusion that since 
all sensors appeared to have operated in an optimum manner, 
and salinity and temperature signals were strong and yielded 
good results, the sound velocity signal was suspect. The 
signal had been filtered and amplified and conceivably could 
have been altered, or else a mistake made in the data re- 
duction settings. In order to determine which values of the 
variance of sound velocity should be used for the remainder 
of the analysis a final check of the results of both methods 
was conducted. Part of the experiment involved the measure= 


ment of sound amplitude modulation (see NPS thesis of LCDR 
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W. J. Smith, Jr., December 1971). Stone and [Ref. 
10] define the coefficient of variation (CV) as: 
Ae 


HE 





CV = 


where AP is the variation of the acoustic pressure relative 
to the mean value P. In his experiment LCDR Smith measured 


Be which he shows is related to CV by: 


(ey) = 40 a,’ 


Chernov [Ref. 11] formulated the relationship: 


few) - = VT po K,- als 
where, 
2T 
Rn = — 
A 


À = acoustical wave length 


= variance of sound velocity 


a = correlation length (Gaussian correlation function 
assumed) 
L = distance between transducer and receiver 


This relationship combined with the measured values of CV 


and sound velocity variance, allowed for the calculation of 


Fie correlation length: 


RUN MT a (from var Ctheory! a (from var Cexpt! 
m) (cm) (cm) 
5 6.9 175 0 2.699 x 10° 
3 9.3 SS = 
4 13.9 TENIS 4279.9 
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The values of "a" calculated from velocimeter data are un- 
reasonably large. For this reason the theoretical sound 
velocity variances are used throughout the remainder of the 
analysis. A graphical plot of theoretical sound velocity 
variances versus time of day with depths in parentheses are 
own in Fig. 15. 

It is interesting to observe that Skudrzyk [Ref. 12] 
also found patch sizes, which increased approximately 
linearly with depth in approximatley the same manner as 
shown above. However, Skudrzyk's patch sizes were much 


larger than reported here, probably due to his acoustic 


range being much greater than in this experiment. 


DAL 





VIE. CONCLUSTEONS 


A. Although the near surface regime is one of complex and 

high variability, several conclusions can be made: 
1. Temperature fluctuations decrease with depth during 
the period just prior to sunset. 
2. Temperature fluctuations increase with depth just 
after sunrise, but are smaller in magnitude than the 
fluctuations just prior to sunset. 
3. Salinity fluctuations decrease with depth but are 
larger in magnitude in the period after sunrise. 
4. Sound velocity fluctuations normally follow the 
trend of temperature variations. However, the im- 
portance of salinity fluctuations cannot be neglected. 
In one instance the salinity fluctuations had as much 
an effect on sound velocity fluctuations as did tem- 
perature fluctuations. 
5. The maximum energy spectra of salinity and temper- 
ature fluctuations occurred at frequencies that agreed 
with dominant surface wave frequencies. 
6. The sizes of patches of sound velocity fluctuations 
were calculated from the statistics of sound amplitude, 
temperature and salrnrtyesluetustrous Tihe implied 
patch sizes increased approximately linearly with in- 


creasing depth. 


SE 





VIII. RECOMMENDATIONS FOR'FURTHER STUDY 


A. While the medium was described in terms of temperature 
and salinity, a major question remains, "Are there other 
variables that could affect the propagation of sound in 
sea water?" In order to best achieve an answer to this 
question one should ensure that a correct experimental 
value of sound velocity is measured for comparison with 
una velocity computed theoretically. LCDR Smith's and 
Rautmann's thesis work indicated the presence of bubbles 
in the medium that could have a marked affect on sound 


velocity. 


B. In order to complete the description of the medium, an 
analysis of the size of the "patches" of temperature and 
salinity should be performed by use of arrays of sensors. 
Application of knowledge of the drift velocity, in addition 
to the procedures used in the present experiment, would 


also enable patch size determination. 
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APPENDIX A 
GENERAL DESCRIPTION OF BISSETT-BERMAN MODEL 9006 SALINITY, 


TEMPERATURE AND DEPTH MEASURING SYSTEM [Ref. 13] 


The model 9006 STD is comprised of underwater sensors 
and deck equipment whose functions are to furnish power to 
the sensors and process data received from same. 

In situ measurement of salinity is determined by measur- 
ing conductivity, temperature and pressure. An inductively- 
coupled sensor enables detection of conductivity, which is 
also compensated for temperature and pressure effects, 
thereby producing an output totally dependent upon salinity. 
The output shifts the frequency of a PARALOC signal pro- 
viding an FM analog of salinity. 

The temperature sensor uses a platinum resistance 
thermometer and a PARALOC which provides an FM analog of 
bemperature. 

The depth sensor utilized a pressure transducer into 
which is incorporated a strain-gauge bridge circuit. The 
changes in resistance are converted to a frequency analog 
in the PARALOC. 

The signal mixer receives and — s power from the 
deck equipment and sends it onto the underwater sensors. 

It then multiplexes and amplified the FM signals from the 
sensor and sends them to the deck equipment. 


The distribution amplifier amplifies the FM signal 


from the underwater units, and separates data signals from 
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the DC power which enters the distribution amplifier from 
the power supply rack. 

Each underwater sensor has a separate discriminator 
which separates the mixed data signal into an individual 
channel for recording availability. The discriminator 
filters the desired signal and rejects others. The accepted 
signal is anplified, squared, and converted to a 10 millivolt 
DC level which is the output of the recorder. Dial readouts 
enable continual monitoring of this signal. 

The Power Supply is a model 8600 and provides a constant 
voltage of 26.5 volt DC to the deck units, and a constant 
current of 150 milliamps at a maximum of 140 volts DC to the 
underwater units. The underwater rack is stainless steel. 
The sensors and mixer are mounted in individual shock- 
absorbent brackets. 

A single conductor cable provides for physical connec- 


tion between sensors and deck equipment. 
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APPENDIX B 


DATA CUMULATION AND RUN DESCRIPTION 


The output DC voltage from the oceanographic sensors 
was recorded on the Sangamo magnetic tape recorder as 
explained in "Experimental Configuration" section. The 


Tis during which data were collected are described below: 


RUN DEPTE DATE START TIME STOP TIME 
1 4. 3m ZU Oct 1419 1429 
2 4.2m 2r OCE 1530 1550 
3 Sum AO 1616 1636 
4 roe 2l OCE 1648 1708 
5 679m AO CE 1728 1748 
6 4. 3m 22 OCE 0354 0448 
7 4.3m ZO CE 0546 0622 
8 Dom DIO EE 0650 0715 
9 14.6m 2208 ct. 0725 0749 

10 "Bem Ż2wWOCE 0802 0821 

11 822m LOSE 0832 0852 

12 5.7m 22 Oct 0929 0945 

SENSITIVITY SETTINGS 
Run 3 
salinity 50 millivolts/millimeter 


thermistor 1 100 millivolts/millimeter 
thermistor 2 100 millivolts/millimeter 


Run 4 
salinity >50 millivolts/millimeter 


thermistor 2 100 millivolts/millimeter 
sound velocity 50 millivolts/millimeter 
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Run 5 


salinity 50 millivolts/millimeter 

thermistor 1 100 millivolts/millimeter 

thermistor 2 50 millivolts/millimeter 

sound velocity 20 millivolts/millimeter 
Run 6 

salinity 50 millivolts/millimeter 

thermistor 1 100 millivolts/millimeter 


sound velocity 50 millivolts/millimeter 


Run 8 
salinity 100 millivolts/millimeter 
thermistor 1 100 millivolts/millimeter 


sound velocity 100 millivolts/millimeter 


Run 9 
salinity 100 millivolts/millimeter 
thermistor 1 100 millivolts/millimeter 
Run 10 
salinity 50 millivolts/millimeter 
thermistor 1 100 millivolts/millimeter 


In addition to Bissett-Berman salinity measurements, 
Nansen casts were taken during each run in order to monitor 
the calibration correction which had to be applied to the 
dial readout of the base value of salinity (see Table III). 
The correction factor was (-) 5.1 ppt. While this was a 
high correction factor, pre-experiment checks on the STD 
led to the assumption that the system was operating correctly 
in the measurement of fluctuations and was consequently used 


for chat purpose. 
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Thermistor 2 Run 3 
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Power Spectral Density 
Salinity Run 3 
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FIGURE 41. 


Power Spectral Density 
Thermistor 2 Run 4 
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Spectrum ((°C)*/HZ) 
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FIGURE 42. 
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Spectrum ((°/,,)°/HZ) 


Versus Frequency (0-2.5 HZ). 
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Power Spectral Density 
Sound Velocity Run 4 
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Spectrum ((m/sec) */HZ) 
Versus Frequency (0-2.5 HZ). 
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Power Spectral Density 
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Spectrum ( (°C) 2 /HZ) 
Versus Frequency (0-2.5 HZ). 
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FIGURE 46. 


Power Spectral Density 
Salinity Run 5 
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Versus Frequency (0-2.5 HZ). 
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FIGURE 48. 


Power Spectral Density 
Thermistor 1 Run 6 

42 3m 0354229606 — 
Spectrum ((?C)^/HZ) 

Versus Frequency (0-2.5 HZ) 
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FIGURE 49. 


Power Spectral Density 
Salinity Run 6 
4.3moO0S5SA2ZZNOGE 

Spectrum 60700752) 

Versus Frequency (0-2.5 HZ). 
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Power Spectral Density. 
Sound Velocity Run 6 
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Spectrum ((m/sec)^/HZ) 
Versus Frequency (0-2.5 HZ). 
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FIGURE 52. 


Power Spectral Density 
Sound Velocity Run 8 

9.5m 0650 22 Oct 

Spectrum ( (m/sec) */HZ) 
Versus Frequency (0-2.5 HZ). 
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Power Spectral Density 
Thermistor Run 10 

7.-8m 0802 22 OCE 

Spectrum ((°C)?/HZ) 

Versus Frequency (0-2.5 HZ). 


| | 





95 








” 
mi 


.. 


dE 


` ewe 
`Q aa 


FIGURE 55. 


Power Spectral Density 
Salinity Run 10 

7.8m 0802322 Oc 
Spectrum ((?/,,)^/HZ) 
Versus Frequency (0-2.5 H2). 
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